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Abstract 

The Young's modulus of a two-dimensional truss-like structure is simulated by us- 
ing the finite element method. A power-law expression is proposed for the effective 
Young's modulus of the system. The obtained numerical results are compared with 
the experimental data of the anisotropic thin cellular polypropylene films. At high 
solid volume fractions (> 0.4), the average shape of the cells are lateral, and their 
dimensions have around one-to-five ratio. As the samples are inflated further, vol- 
ume fraction of the solid is decreased, the average shape approach a diamond-like 
structure with one-to-two ratio. In addition the effective Young's modulus of the 
system increases. It is concluded that valuable structural information can be ob- 
tained by analyzing the experimental data and the numerical simulations, which 
take into account the material's micro-structural information, simultaneously. 
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1 Introduction 



After the invention of the first man-made cellular material in 1930's by Munters 
and Tandberg |l| , these new high-tech materials were introduced in insulation 
and packaging products in the daily life. Nowadays, these lightweight ma- 
terials, in the form of voided-solids or foams, are widely used in variety of 
engineering applications [2, || such as thermal insulation, aeronautics, impact 
absorbing materials g [J dectro-mechanical sensors |6(, etc. In this paper, 
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we present the mechanical properties (effective Young's modulus) of a layered 
cellular structure, in which the layers form a truss- like micro-structure. The 
macroscopic properties of the model structure are calculated using the finite 
element (fe) method 0, H, E| EH- The effective Young's modulus of the 
considered structure can be expressed by a power-law function similar to the 
one proposed by Gibson and Ashby j3| • Here, the expont in the power-law is ex- 
pressed as a function of solid volume fraction q. The computed values are com- 
pared with the experimental results of thin cellular polypropylene films ^ . 



It is shown that there is a good agreement between numerical calculations and 
experimental measurements. The comparison also gives valuable information 
regarding micro-structure of the polypropylene films. The information could 
be used to optimize the materials' property for specific applications. The pa- 
per is divided into sections. In the proceeding section (§ 2), a brief review on 
homogenization is explained. In § 3 the geometrical and computational models 
are presented. In addition, the numerical results, and the influence of various 
geometrical parameters are presented and discussed in the same section. The 
numerical results are compared to those of experiments in § 4. Conclusions 
are in 8 5. 



2 Background 

The foam structure is a heterogeneous medium whose macroscopic or effective 
properties can be expressed as a function of its constituents. 

where *p e is the effective property sought and subscript 'i' represents the prop- 
erties the constituents with volume fractions q^ YJ Qi — 1 (in Gibson and Ashby 
0] notation q = p* / p, where p and p* are the density of the solid and foam, 
respectively). Determining *}3 e for heterogeneous materials or structures or 



finding a representative function / is an old problem in science [13J, |14| . The 



simplest forms for function / are the arithmetic f a and harmonic fh averages, 

/.(fc,*Pi) = Efc*Pi and ( 2 ) 
A(ft,Vi) = (Eft/Vi) _1 . (3) 

These two equations are actually the bounds for mixture properties [l5l . fl6L fl7j ] . 
The former one, Eq. (2), represents a mixture with layers parallel to the applied 
field. The same arrangement perpendicular to the applied field direction is 
expressed by Eq. (3). For complex mixtures, i.e. the constituents are arranged 
in a randomly fashion, non-trivial mixture formulas and rigorous bounds are 
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applied [HQ. Eq. (2) is often valid for low volume fractions of constituents, 
however, for open-cell cellular solids, the effective Young's modulus E e is not 
trivial to find due to lack of mechanical properties of the void. It was stated 
that E e is related to the volume fraction of the solid material q through a 
power-law relation 0], 



E e = £q n E s (4) 

where, £ and n depend on the micro- structure, and E s is the Young's modulus 
of the solid. The value of n lies in the range 1 < n < 4, yielding a wide range of 
effective properties for a given solid volume fraction q 19]. Gibson and Ashby 
0] have presented experimental results on various cellular materials, and their 
observations have suggested n = 2. Since the dependence of effective properties 
on micro-structure are not well understood in mixtures, the exact form of £ 
and n are not known, which is an handicap to optimize and predict com pos ite 
properties. However, computer simulation are prefered in such cases 



>mposi 



3 Model cellular structure 



3.1 Geometrical consideration 



The considered two-dimensional truss-like structure is presented in Fig. 1. It 
is a symmetric geometry with solid regions and diamond-like voids centered 
at the triangular lattice nodes. Because of the lattice structure, the region 
marked between [0, a/2] and [0, 6/2] in Fig. 1 is used as the unit cell in 
the simulations. The marked region is the smallest representative part of the 
whole structure with a/2 and 6/2 in the x- and ^-directions, respectively. The 
geometrical model of the unit cell is illustrated in Fig. 2. The dark polygon 
between ABCDEF represents the solid medium Q s with wall thickness t. The 
expressions that define the geometrical parameters are as follows. 



q = ah /4 — mn, (5) 
a = arctan(6/a), (6) 



m = y (1 — g)a6/(4tana), (7) 
n = m tan a, (8) 

t(q) = 2^(6/2 - n) 2 + (a/2 - m) 2 sin(2a) (9) 

where q is the area/volume fraction of the solid medium and = 2a in Fig. 2. 
The angles a and are defined as illustrated in Fig. 2, and are used to calculate 
the cell-wall thickness t. The introduced parameters m and n are related to 
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the dimensions of the diamond-like void. From the above equations, one can 
also calculate the concentration q for a given thickness t and the unit-cell 
dimensions a and b, 



<^- a( ' + ^; + -y- ( > T) • ■ m> 

where r = t csc(2a). In the FE analysis, the effective Young's moduli are 
calculated for various a and b values. 



3.2 Numerical analysis 



Our FE method is based on the minimization of the elastic energy in con- 
junction with adaptive meshing 0,H^. The Newton's second law in the two- 
dimensions, which is in reality a system of two linear equations known as 
Navier's equations, are solved for spatial points obtained from the FE trian- 
gulation of Q s . 



-V T = K (11) 

Here T is the mechanical stress tensor and K is the load vector. The stress 
tensor has three components in two-dimensions, which are normal stress com- 
ponents T x and T^, and the shear stress component S^, 



ox ay 

-nf-t^ < 13 > 

Including the Hooke's law between the stress and the elastic strain, T = c Vu, 
lead to a partial differential equation expressed in the global displacements 
u(u, v). Then, Eq. (11) becomes as follows 



- V ■ cVu = K, (14) 

where the coefficient c is a function of solid materials Young's modulus and 
Poisson's ratio. In Fig. 2, the shaded areas and the arrow represent the con- 
straints and load, respectively. The Young's modulus of the structure is calcu- 
lated from the average deformation of the edge AB in the ^-direction / vdx/(l— 
m). 
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(15) 



The solid region fi s in Fig. 2 is used in the simulations together with the 
boundary conditions presented in Table 2. Meshing of the computational do- 
main is an important factor in the FE analysis, we have, therefore, adapted 
a procedure that uses the cell-wall thickness as the meshing parameter. In 
the calculations, the FE triangle sides are not larger than the one-tenth of the 
cell-wall thickness. This procedure allowed us to discretize even very thin cell 
walls. 

In the calculations the Young's modulus of the solid is E s = 1 GPa 21], and 
its Poisson's ratio is v s « 0. The effective Young's moduli of the structures are 
practically independent of the u s , which has also been observed previously [HI 



3.3 Illustrative examples 

Three simulations are illustrated in Fig. 3. The gray scale map and contours 
represent the global deformation and von Mises stress distributions T v m, 



In the figures, the dimensions of the representative units are varied in order to 
display the change in the stress and displacement distributions for selected a/b 
and q values. For given structural parameters a/b and the volume fraction of 
the solid q yield distint global displacements and von Mises stress. The effective 
Young's moduli for these structures are 170 MPa, 3.4 MPa and 1.0 MPa for 
Figs. 3a, 3b and 3c, respectively. The positions where the max and the min 
deformations are observed are also indicated in the figures. Depending on a/b 
and q, the bending of the beam-like cell-wall structure is peculiar, and the 
cell-wall thickness is an important factor for this behavior. 

3.4 Effective properties of the cellular structure 

The effective Young's moduli, calculated from the deformations, for various 
a/b ratios are shown as log-log and semi-log plots in Fig. 4. In the figure, four 
of the simulation results are illustrated, a/b = {1, 2, 4, 8}. At low concen- 
trations of the solid (q < 0.20), the effective Young's modulus E e is a linear 
function of the volume fraction in the log- log representation [as in Eq. (4)]. 



23. 




(16) 



5 



Table 1. 



Fitting parameters of Eq. (17) at different volume fraction intervals, the sectors displayed in Fig. 4. 





Q 


< 0.20 




0.20 < 


: q < 0.55 


0.55 < q < 0. 


80 




q > 0.80 




a/b 


log Oi 


bi 


Ci 


log a 2 


b 2 


c 2 






c 3 




b 4 


c 4 


1 


-0.69 


-0.67 


0.96 


-0.14 


1.53 


1.30 


0.93 


2.76 


1.04 


0.99 


4.61 


-0.14 


2 


-1.50 


-0.77 


0.96 


-0.46 


3.59 


1.55 


1.08 


8.96 


0.49 


1.00 


4.66 


3.59 


3 


-2.10 


-0.77 


0.96 


-0.82 


4.70 


1.68 


1.13 


14.21 


-0.21 


1.01 


5.81 


6.02 


4 


-2.56 


-0.79 


0.96 


-1.15 


5.32 


1.74 


0.95 


17.62 - 


-0.71 


1.03 


12.91 


3.45 


5 


-2.93 


-0.77 


0.96 


-1.45 


5.63 


1.77 


0.70 


19.63 - 


-0.97 


1.05 


24.65 


-3.12 


6 


-3.24 


-0.80 


0.96 


-1.71 


5.87 


1.79 


0.49 


20.85 - 


-1.11 


1.06 


38.87 - 


-11.99 


7 


-3.49 


-0.77 


0.96 


-1.94 


6.01 


1.81 


0.34 


21.64 - 


-1.19 


1.06 


53.98 - 


-21.87 


8 


-3.72 


-0.79 


0.96 


-2.14 


6.12 


1.82 


0.23 


22.13 


-1.22 


1.05 


68.86 - 


-31.84 


9 


-3.93 


-0.84 


0.95 


-2.32 


6.22 


1.82 


0.17 


22.51 - 


-1.24 


1.02 


83.45 - 


-41.81 


10 


-4.11 


-0.84 


0.95 


-2.50 


6.20 


1.83 


0.12 


22.70 - 


-1.21 


0.99 


97.34 - 


-51.42 



However, the constant <£ is related to the unit-cell dimensions a and b. When 
the volume fraction of the solid is increased there is no simple relation between 
E e and q, which is illustrated in Figs. 4a and 4b. In the intermediate low con- 
centrations 0.20 < 0.55, the dependence of E e is slightly semi-logarithmic 
with similar slopes but the shifts are different related to the unit-cell dimen- 
sions (Fig. 4b). At higher concentration of solid, the behavior of E e is still a 
power-law dependence on q. It is non-trivial to find an homogenization for- 
mula that can express the effective Young's modulus E e as a function of whole 
range of solid volume fraction (0 < q < 1) and the unit-cell dimensions for the 
truss-like structure. However, by dividing the g-axis into four separate sectors. 
The facing sectors are separated at q = {0.20, 0.55, 0.80}, and the dashed ver- 
tical lines (----) in Fig. 4 are drawn to represent the borders. We are able 
to represent the effective Young's modulus of the truss-like structure with a 
power-law similar to Eq. (4) in each sector when we modify the exponent u(g), 
n = aq + b. The proposed homogenization equation is then expressed as 



Et(q) = [a iq ^ q+u) ]E s with 2 = 1,2,3,4 (17) 

where dj, bj and Cj are the fitting parameters, and i is the sector number. The 
values of these parameters obtained from a curve-fitting procedure for each 
sector are presented in Table 1. The resulting curves are shown in Fig. 4 with 
solid lines (— — ). It is remarkable that only CI4, b\ and Ci show very slight 
dependence-they are nearly constant-on the dimensions of the representative 
unit a and b, and bi q + Ci ~ 1. Last but not least, the curves are shifted down 
with increasing a/b, because of the increasing anisotropy in the structure. 

The behavior of effective properties are of importance at low concentrations 
of solid material for technological applications, since the material would be 
light (the density of solid would be low) For this reason, we first pay 

our attention to sector 1. In this sector, the exponent term in Eq (17) is 
lower than 2 and it is approximately 1, which coincides with the arithmetic 
averages, Eq. (2). Moreover, when n ~ 1 , it is an indication that the cell- 
wall stretching is a dominant behavior 0, IHI]. Here, it is shown that for two 
dimensions this exponent is also 1, which clearly illustrates that we simulate 
the cell-strecthing behavior for q < 0.2. Although, E e /E s oc q, one should 
adopt an effective volume fraction, di in our case, as a function of structure 
dimensions a and b. A sum of an exponential and a power-law was sufficient 
to describe the behavior of Oi, 



Oi(a, b) = A exp[B (a/b)} + C (a/b) B (18) 

Here the constants A = —16.84 B = —3.845 and C = 0.565, respectively. As 
the ratio of the cell dimension increases a/b — > 00, Oi(a/6) — > 0. As a result, 
for low concentrations of solid, the effective Young's modulus of the truss-like 
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structure in two- dimensions can be expressed as 

Ef(q < 0.2) « {A exp[5 (a/b)] + C {a/b) B } q (19) 



In sectors 2, 3 and 4, the relations between fitting parameters and a/b are non- 
trivial. Although the fitting parameters in Table 1 indicate clear dependence 
on the unit-cell dimensions a/b, expressions in the form of Eq. (18) are not 
sought. 



As a note, if we consider the similarity to the Bruggeman symmetrical |24L l2 
formula or general Wiener formula j2y, 27 , 28| , which use a dimensionality pa- 
rameter d in the effective permittivity expressions, the pre-exponent cti here is 
related to the internal structure of the material at low concentrations of solid. 
Tuncer et al. 26|| have stated that the shape parameter in dielectric mixtures 



could in principle be used as the dimensionality. In the present problem, the 
shape parameter is hidden in the a/b, such that when a/b — > oo the structure 
is effectively a lateral layered structure, and as stated in § 1, the void does not 
have a mechanical property (E v = 0), therefore E e — > even for high vol- 
ume fractions of solid medium. At high a/b values (for anisotropic structures), 
the effective load carrying volume of the solid is small. This is also observed 
in Fig. 4 in sectors 3 and 4 on which E e /E s drops fast with increasing a/b. 
And when a/b — ► 1, the change in E e is not as drastic as a/b > 1 cases as 
q decreases from 1 to 0.6, due to the effective load carrying sections of the 
structure are larger in volume fraction. 



3. 5 Influence of the cell-wall thickness 



Gibson and Ashby [2j has shown that the mechanical properties of foams could 
be scaled when the cell-wall thickness was much smaller than the relative cell 
size (q 1), 



E e = <B(t/b) m E s , (20) 

where 03 and m are scaling parameters, and t and b are the cell-wall thickness 
and unit-cell height for a/b = 1. They have stated that the value of m is 3 and 
4 for open- and closed cell foams, respectively. 

The effective Young's modulus as a function of the cell-wall thickness is illus- 
trated in Fig. 5 for 10 different unit-cell dimensions. The behavior is obtained 
by converting the relation in Eq. (17) as a function of cell-wall thickness by 
substituting Eq. (10) in Eq. (17). It is clearly indicated that the thicker walls 
yield greater mechanical stiffness because of the higher volume fraction of the 
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solid q [29|. This result is also discussed in § 3.4. Structures with thin walls 
and low q solid yield very low values of elastic moduli, and E e /E s oc (t/a). 

There is a critical knee-point t c at which the power-law behavior of E e changes 
its exponent from a lower value to a higher one. As we approach the criti- 
cal thickness from high thickness values the physical nature of the truss-like 
structure changes the stress field distribution, and the deformation is altered. 
The position of the critical thickness is calculated from a similar expression 
as in Eq. (17), however, only two regions are considered, t/a < t c /a and 
t/a > t c /a-the regions are separated by t c /a. The logarithm of the critical 
thickness against (a/b) 1 / 2 indicate linear dependence on the unit cell dimen- 
sions, which is shown in the inset (a) of Fig. 5 and E e /E s oc (t/a) ' 
If a similar relation as Eq. (20) is considered for t/a > t c /a, the exponent in 
Eq. (20) has an a/b dependence, m = (a/b). Finally, since a/b is related to the 
angle between the truss-like surfaces, for a given volume fraction q, the soft 
and stiff structures can be tailored by varying the unit-cell dimensions. This is 
illustrated in the inset (b) in Fig. 5. At large angle values a > 20°, the critical 
thickness changes slowly compared to the values at low angles a < 15°. 



4 Comparison of numerical results with experiments 



In Fig. 6, the scanning electron microscope cross-sectional view of a cellular 
polypropylene film is shown. The dark areas indicate the void regions, and the 
bright regions are the solid polypropylene, Q s . The original cellular polypropy- 
lene films were supplied from VTT Processes, Finland. They were 37 (xm thick 
and stretched [the detailed description of the material can be found in Ref. 
3oj ] . The films were inflated bv means of a pressure treatment in order to vary 
the dimensions of the cells The treatment not only influences 



the shape and the size of the cells, but also the relative density (the solid 
volume fraction q). Several samples were prepared and inflated for each indi- 
vidual measurement. The effective Youn g's modulus of the samples has been 
estimated from dielectric measurements 34 



The micro-structure of the polypropylene films indicate a cell structure similar 
to the regular geometry used in the FE analysis in § 3.1, but the cell- walls 
do not show a sharp connection points (corners) as in Fig. 1 (or the nodes 
E and F, and B and C in Fig. 2). The structure show high anisotropy, and 
the cell-dimensions are not equal, a/b ^ 1. Moreover, the size distribution 
of voids in the polypropylene are not constant but distributed-the structure 
is irregular. Yet, comparison with the simple truss-like model would yield 
the most probable or in other words the average void size and the cell wall 
thickness for material design purposes. 
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To display the utility of the FE modeling, the simulations are compared with 
experimental data of the cellular polypropylene films. The numerical data and 
the experimental results are ploted against the volume fraction of the solid 
(polypropylene) q in Fig. 7. Iso-thickness lines from the simulations are also 
displayed in the figure. The experimental data are shown with symbols (□) 
and error bars (the samples were thin (< 100 yim) and soft, therefore, the 
uncertainty in the thickness determination leads the error bars). The Young's 
moduli of polypropylene samples were normalized by the value for the extru- 



sion grade polypropylene [21|, E s w 1 GPa. The measured values were between 
1 MPa and 10 MPa depending on the density of the samples. The logarithm 
of the data indicate a quadratic dependence on the volume fraction 



\ogE e /E s = Wx 2 -8.9x-0.84. (21) 

The dashed line (----) corresponds this expression in the figure. In order 
to show the agreement between the numerical and experimental behavior, we 
apply the same curve fitting procedure to the iso-thickness line at t/a = 0.075. 
Similar to the experimental data, the considered iso-thickness line can be 

expressed by a quadratic expression, shown as the dotted line ( ) in the 

figure. This illustrates that the change in the behavior of stiffness might be 
due to the alteration of the average cell dimensions. 



However, a detailed analysis of the experimental data and the numerical sim- 
ulations indicates that by inflating the samples first the cell-wall thickness is 
altered-cell-wall stretching 0,H3]. The experimental data on q > 0.4 can be ex- 
pressed by the truss-like structure with cell dimensions a/b w 5. When q < 0.4, 
the data bend, and the sign of dE e /dq becomes negative. The Young's moduli 
of the samples slightly follow the iso-thickness line t/a = 0.05 for q < 0.4. 
This shows that the cell-wall stretching is decreased. 



The critical thickness for polypropylene is around 0.05. The samples become 
stiffer at low solid concentrations. The extrapolation of the data using Eq. (21) 
yields the thinnest cell-wall thickness t/a pa 0.03 for the polypropylene films. 



In overall, the experimental data lie inside the region marked by2<a/6<6 
and normalized thickness 0.03 <t/a< 0.12. Since, monitoring and measuring 
the Youngs'moduli of the samples in-situ are not presently possible during 
inflation, we have not been able to indicate the experimental traces of the 
iso-thicknesses. However, the results indicate that there is a good aggrement 
between the behavior of the inflated samples and the model structure used in 
the FE analysis. 
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5 Conclusions 



In this paper, we present the numerical calculation of effective Young's modu- 
lus for a two-dimensional truss-like structure. A new homogenezation formula 
is proposed, which expresses the effective properties of the system in a broder 
concentration range. In this procedure the concentration axis is divided into re- 
gions. The expression yields a log-log dependence of effective Young's modulus 
on solid concentration q when g_< 0.2. Similar kind of behavior has also been 
observed for foams previously 0, Hif . The numerical methods were compared 
with experimental data of thin cellular polypropylene film. The comparison 
shows that once the density of the inflated samples are known, we are able to 
estimate the average cell dimension and cell-wall thickness. The samples in- 
dicated that cell-wall stretching is important at solid concentrations over 0.4, 
which yield cell dimensions a/6 ~ 5. At high solid concentration, most of the 
solid material are mechanically passive yielding low effective Young's modulus. 
As the samples were inflated even further (i) the cell-dimensions are altered, 
(ii) stiffer structures with a/b ~ 2 are obtained, and (Hi) the cell-stretching 
behavior is decreased. As a final remark, the manufacturing process of the 
cellular polypropylene films picks out particular cell-wall thicknesses. 
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Fig. 1. Considered truss-like geometry in the cellular material simulations. The 
shaded area represents the smallest representative unit (unit-cell). 
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Fig. 2. Unit cell (representative volume element) in the simulations. The dark-gray 
region (ABCDEF polygon) is the solid medium, Q s . The shaded areas represent 
the boundary conditions, such that the solid region is allowed to deform in the 
rectangular area marked by AHDG. The arrow represents the force (load) applied 
to the edge AB. The wall thickness of the solid is calculated using the angles GEF 
and EFB, which are denoted by a and 4>, respectively. The volume fraction of the 
solid region Q s , polygon ABCDEF, is ab/4 — mn. 
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Fig. 3. Global displacement (gray scale color map) and von Mises stress (contour 
plot) distributions for three selected structures (a) a/b = 1 and q = 0.5, (b) a/b = 4 
and q = 0.5 and (c) a/b = 4 and q = 0.25, respectively. Darker regions show the 
highest deformations. The minimum and maximum deformation locations for unit 
applied loads are indicated in the graphs in metric units. The calculated effective 
Young's moduli for these structures are (a) 170 MPa, (b) 3.4 MPa and (c) 1.0 MPa, 
respectively. 
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Fig. 4. Normalized effective Young's modulus E e /E s as a function of volume fraction 
of solid material q, (a) in log-log and (b) in semi-log plots. Only four a/b ratios are 
shown. The curves are divided into four sectors as separated by dashed vertical lines 
(----) at q = {0.2, 0.55, 0.8}. The effective Young's modulus in each sectors are 
expressed with a power-law, Eq. (17). The solid lines (— — ) represent the power-law 
curves with the fitting parameters presented in Table 1. 
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Fig. 5. Normalized effective Young's modulus E e /E s as a function of normalized cell 
wall thickness t/a. The critical thickness for a/b = 10 is marked with t c , and the 
inset (a) displays the linear dependence of log(i c /a) on square root of dimensions 

ratio (a/6) 1 / 2 , the solid line ( ) is log(t c /a) = 0.32 (a/b) 1 / 2 - 0.4. The inset (b) 

is the same dependence as a function of the angle a defined in Eq. (6) 
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Fig. 6. Cross-sectional scanning electron microscope image of a cellular polypropy- 
lene film. The dark regions are air-voids and bright spots are polymer. The film 
was inflated with 25 bar at 90 °C from a 0.73 g/cm 3 and 37 |J.m thick sample. The 
density of the sample is 0.44 g/cm 3 . 
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Fig. 7. Normalized effective Young's modulus E e /E s as a function of solid volume 
fraction q. The experimental data are illusrated with symbols (□). Iso-thickness 
lines are also shown with thick lines, and corresponding normalized thickness values 

are indicated for each line. The dashed (----) and dotted ( ) lines represent 

quadratic relations for the experimental (log E e / E s = 10 x 2 — 8.9 x — 0.84) and the 
iso-thickness line t/a = 0.075 (log E e / E s = 9.4x 2 — 11 x + 0.59), respectively. 
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Table 2 

Boundary conditions in the FE. The constrained edges have u = or v = conditions 
for 'yes' labels. All the other edges are 'free' to move in corresponding directions. 
The a pplied load per unit length is presented by Fp. 

AB BC CD DE EF FA 

Load F 

x-constraint yes free yes yes free yes 

y-constraint free free free yes free free 
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